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Abstract 

We consider the linear stochastic wave equation with spatially homoge- 
nous Gaussian noise, which is fractional in time with index H > 1/2. 
We show that the necessary and sufficient condition for the existence of 
the solution is a relaxation of the condition obtained in [10], when the 
noise is white in time. Under this condition, we show that the solution 
is L 2 (f2)-continuous. Similar results are obtained for the heat equation. 
Unlike the white noise case, the necessary and sufficient condition for the 
existence of the solution in the case of the heat equation is different (and 
more general) than the one obtained for the wave equation. 

MSG 2000 subject classification: Primary 60H15; secondary 60H05 
Keywords and phrases: stochastic wave equation, random field solution, spa- 
tially homogenous Gaussian noise, fractional Brownian motion 

1 Introduction 

The random field approach to s.p.d.e.'s initiated in [IB], has become increasingly 
popular in the past few decades, as an alternative to the semigroup approach 
developed in [T^ , or the analytic approach of [5U] . 

Generally speaking, a random field solution of the (non-linear) equation: 

Lu(t,x)=a(u(t,x))W(t,x)+(3(u(t,x)), t>0,xeM. d (1) 
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(with vanishing initial conditions) is a collection {u(t,x),t > 0,x G M. d } of 
square integrable random variables, which satisfy the following integral equa- 
tion: 

u(t,x) = / G(t-s,x-y)a(u(s,y))W(ds,dy)+ / G{t-s,x-y)(3(u(s,y))dyds 

JO JR d JO JR d 

provided that both integrals above are well-defined (the first being a stochastic 
integral). In this context, L is a second-order partial differential operator with 
constant coefficients, G is the fundamental solution of Lu — 0, and W is a 
formal way of denoting the random noise perturbing the equation. 

When the equation is driven by a space-time white noise (i.e. a Gaussian 
noise which has the covariance structure of a Brownian motion in space-time), 
the random field solution exists only if the spatial dimension is d = 1. In 
this case, the stochastic integral above is defined with respect to a martingale 
measure, and the solution is well-understood for most operators L, in particular 
for the heat and wave operators (see [46], [7J or [43]). 

To obtain a random field solution in higher dimensions, one needs to consider 
a different type of noise, which can be either Gaussian, but with a spatially 
homogenous covariance structure given formally by: 

E[W(t,x)W(s,y)]=S(t-s)f(x-y), 

or of Poisson type. Historically, the two approaches have been initiated at about 
the same time (see [27], [TT], [25] for the wave equation with Gaussian noise in 
dimension d = 2, and [5], [U for the Poisson case). 

After the ingenious extension of the martingale measure stochastic integral 
due to |10| . it became clear that the random field approach can be pursued for 
the study of general s.p.d.e.'s with spatially homogenous Gaussian noise. Since 
this extension allows for integrands which are non-negative measures (in space), 
the theory developed in [10] covers instantly the case of the (non-linear) wave 
equation in dimensions d E {1,2,3}, and the case of the heat equation in any 
dimensions d. In the non-linear case, the existence of the solution is obtained 
by a Picard's iteration scheme, under the usual Lipschitz assumptions on a, f3, 
and the following condition, linking the operator L and the spatial covariance 
function /: 

/ / \TG(u,-m\ 2 d Uf i(dO<™. (2) 

JR d JO 

(Here [i is a non-negative tempered measure, whose Fourier transform in /.) 

Moreover, ([2]) is the necessary any sufficient condition for the stochastic in- 
tegral J J Rd G(t — s , x — y)W (ds , dy) to be well-defined, and hence the necessary 
any sufficient condition for the existence of the solution in the linear case, when 
a = 1 and = 0. Since for both heat and wave operators, 

cPt^ < [\TG(u^)(0\ 2 du<ci 2) —^—, forall£eM d , ( 3 ) 
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for some constants , Cj > 0, condition ([2]) is equivalent to: 

—^V(d£) < oo. 

Subsequently, using the Malliavin calculus techniques, it was shown that 
the random variable u(t, x) has an absolutely continuous law with respect to 
the Lebesgue measure on K, and this density is infinitely differentiable. These 
results are valid for the heat equation in any dimension d, and for the wave 
equation in dimension d G {1,2,3} (see |3B], [37], B3J)) under the additional 
assumption (which was removed in [30]): 

L d ( i+W ) < °°' for some a g ^ ^ 

Under ([5]), one also obtains the Holder continuity of the solution for the heat 
equation in any dimension d and the wave equation in dimensions d S {1, 2, 3}. 
This is done using Kolmogorov's criterion and some estimates for the p-th mo- 
ments of the increments of the solution (see [35], [ID], [T2"]V 

The case of the wave equation in dimension d > 4 was solved in the recent 
article [5] , using an extension of the integral developed in [TU] . The existence of 
a random-held solution is obtained under condition @. In the the affine case 
(i.e. a(u) — au + b,a,b G R and /3 = 0), and under the additional assumption 
(J2|), the solution is shown to be Holder continuous. 

In parallel with these developments, a new process began to be used inten- 
sively in stochastic analysis: the fractional Brownian motion (fBm) with index 
H G (0, 1), a zero-mean Gaussian process (-B*)t>o with covariance: 

R H (t, S )= 1 -(t 2H + S 2H -\t- S \ 2H ). 

The case H = 1/2 corresponds to the classical Brownian motion, while the 
cases H > 1/2 and H < 1/2 have many contrasting properties. We refer the 
reader to the survey article [5H] and the monographs [6J and {2S] for more details. 
Most importantly, in the case H > 1/2, 

R H (t,s)=a H [ f \u-v\ 2H - 2 dudv, (5) 
Jo Jo 

where an = H(2H — 1). This shows that {B t )t>o has a homogenous covariance 
structure, similar to the spatial structure of the noise W considered above. 

Returning to our discussion about s.p.d.e.'s with a Gaussian noise, it seems 
natural to consider equation ^ , when the covariance of the noise W is given 
formally by: 

E[W(t,x)W(s,y)] =a H \t-s\ 2H - 2 f(x~y). (6) 

However, this simple modification changes the problem drastically, since 
unless H = 1/2, the fBm is not a semimartingale, and therefore the previous 
method, based on martingale measure stochastic integrals, cannot be applied. 
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Several methods have been proposed for developing a stochastic calculus with 
respect to ffim: (i) the Malliavin calculus (see [14], pQ, [2], [29]), which exploits 
the fact that the fBm is Gaussian; (ii) the method of generalized Lebesgue- 
Stieltjes integration (see [17]), which uses the Holder continuity of the fBm 
trajectories; (iii) the rough path analysis (see [55], [53]), which uses the fact 
that the paths of the fBm have bounded p- variation, for p > 1/H; (iv) the 
stochastic calculus via regularization based also in general on the properties of 
the paths of the fBm (see [TB]V 

These methods have been applied to s.p.d.e.'s (see [53], [3T], [35], [XT]). 
[35]), but not using the random field approach A notable exception is the heat 
equation. The linear equation with noise ^ and H > 1/2 was examined in 
[3J, for particular functions / (e.g. f(x) = with a € (0, d)). We also 

mention the works [32] and [44] for the case of the space variable belonging to 
the unit circle. The quasi-linear equation (i.e. a = 0) was treated in [33], and 
the equation with multiplicative noise (i.e. a(u) = u, (3 = 0) was studied in [15] : 
in these two references, the covariance structure of the noise is a particular case 



(This type of noise is called fractional Brownian field.) The heat equation 
with multiplicative noise © was studied in [1] (for particular functions / and 
H > 1/2) and [19] (in the case H £ (0, 1) and f — So). In the case when the 
spatial dimension is d = 1, the non-linear equation has been treated in |38j using 
a two-parameter Young integral based on the Holder continuity of fBm. 

To the best of our knowledge, there is no study of the wave equation driven 
by a noise W, whose covariance is given by ©■ The goal of the present article 
is to start filling this gap, by identifying the necessary and sufficient conditions 
for the existence of a random field solution of the linear wave equation with 
noise (O and H > 1/2. We also treat the heat equation. 

When H > 1/2, it turns out that under relatively mild assumptions on the 
fundamental solution G of the operator L, the necessary and sufficient condition 
for the existence of the random-field solution of the linear equation Lu = W is: 



which is more general than ([2]). Note that the integrand of the n(d£) integral in 
{7J is the H(0, i)-norm of the function u JFG(w, •)(£)■ Quite surprisingly, and 
in contrast with (J3J , the estimates that we obtain for this norm are different in 
the case of the wave and heat operators: in the case of the wave equation, {7J 
is equivalent to 



of ©: for H,Hi > 1/2 



d 



E[W(t, x)W(s, y)] =a H \t- s\ 2R - 2 ]\{a Hi \x, - yi \ 2H ^ 2 ). 
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whereas in the case of the heat equation, ([7]) is equivalent to: 



) 



2H 



£t(d£) < oo, 



(9) 



The amazing fact is that for the wave operator, these estimates can be de- 
duced using only the estimates of the L 2 (0,t)-norm (given by (J3j) ) , the trick be- 
ing to pass to the spectral representation of the 7i(0, t)-norm of u i— > TG(u : •)(£)■ 
In the case of the heat operator, there is no need for this machinery, since 
u i ► J-G(u, •)(£) is a non-negative function, and its 7i(0, i)-norm can be bounded 
directly by the L 1/ff (0,i) -norm, which is easily computable. 

This article is organized as follows. Section 2 contains some preliminaries, 
and a basic result which ensures that under (|7|) , the stochastic integral of the 
fundamental solution G of the wave operator is well defined. In Section 3, 
we show that the solution of the wave equation exists if and only if ([5]) holds 
(Theorem 13. ip . Moreover, the solution is L 2 (f2)-continuous. Similar results are 
obtained in Section 4 for the heat equation, using (|9|). Appendix A contains 
some useful identities, which are needed in the sequel. Appendix B gives the 
spectral representation of the 7i(0,£)-norm of the function sin. 

2 The Basics 

We denote by C^°(M. d+1 ) the space of infinitely differentiable functions on M. d+1 
with compact support, and <S(R d ) the Schwartz space of rapidly decreasing C°° 
functions in R d . For <p £ L 1 (R d ), we let Tp be the Fourier transform of ip: 



We begin by introducing the framework of [10 . Let fi be a non-negative 
tempered measure on K d , i.e. a non-negative measure which satisfies: 



Since the integrand is non- increasing in I, we may assume that / > 1 is an 
integer. Note that 1 + |£| 2 behaves as a constant around 0, and as |£| 2 at oo, 
and hence (TlO|) is equivalent to: 





/i(<i£) < oo, for some I > 0. 




and 




for some integer I > 1. (10) 



Let / : R d 



R+ be the Fourier transform of /i in <S'(M d 



), i.e. 
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Simple properties of the Fourier transform show that for any tp,ip G 



<p(x)f(x - y)ip(y)dxdy = / ;Fip(£).F^(£)/i(<2£) 



An approximation argument shows that the previous equality also holds for 
indicator functions ip = 1a, ip = Is, with A, B G £>b(R d ), where Bb(M. d ) is the 
class of bounded Borel sets of M. d : 

f(x-y)dxdy= f ^U(e)JMl)M^)- (11) 

JR d 
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As in [3], U, on a complete probability space (fi, J 7 , P), we consider a zero- 
mean Gaussian process W = {Wt(A);t > 0, A G Bb(M. d )} with covariance: 

E(W t (A)W s {B))=R H (t,s) [ [ f{x-y)dxdy=:{l mxA ,l %s]xB ) HV . 

J A J B 



Let £ be the set of linear combinations of elementary functions l[o,t]xA> 
t>0,Ae B b (R d ), and HP be the Hilbert space defined as the closure of £ with 
respect to the inner product (•, -)nv- (Alternatively, HP can be defined as the 
completion of C^°(R d+1 ), with respect to the inner product (•, -)ht-) 

The map l[o,t] xA l— * W t {A) is an isometry between £ and the Gaussian space 
H w of W, which can be extended to HP. We denote this extension by: 

ip^W(ip)=[ [ <p(t,x)W(dt,dx). 

JO JR d 

In the present work, we assume that H > 1/2. Hence, ([5]) holds. From (jlip 
and ([5]), it follows that for any tp,ip G £, 



oo poo 



(<P,iP)hv = olh I J I / ip(u,x)ip(v,y)f(x -y)\u-v\ 2H 2 dxdydudv 



OO POO 



= a H / / Tcp{u,-m^{v,-){i)\u-v\ iH -^{d^)dudv. 

JO JO JR d 

Moreover, we can interchange the order of the integrals dudv and n{d£), 
since for indicator functions ip and ip, the integrand is a product of a function 
of (u, v) and a function of £. Hence, for (p,ip £ £, we have: 

(<p,i>) HV = a H / / ^(u, -XO-^K -)(0W-v\ 2H - 2 dudvn(d0. (12) 

JR d JO JO 

The space HP may contain distributions, but contains the space \HP\ of 
measurable functions (p : R+ x M d — > K such that 

2 / / / / !,„/„, _M|,„/„, „\|i/„ „,\|„, „,|2_ff-2. 



L P\\\HT\ a H I III l^( u ' x ) 1 1^(^12/) 1/(^-2/) h"— "^1 dxdydudv<oo. 

Jo Jo JR d JR d 

We recall now several facts related to the fBm (see e.g. [28]). 
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Let B = (B t ) t >o be a fflm of index H > 1/2. For a fixed T > 0, let H(Q,T) 
be the Hilbert space defined as the closure of £(Q,T) (the set of step functions 
on [0, T]), with respect to the inner product: 

(l[o,t]j l[o,«])w(o,T) = s )- 

One can prove that 

/•Ms 

R H (t,s)= / K H (t,r)K H (s,r)dr, 
Jo 

where A>(t,r) - £(« - r)"" 3 / V" 1 /^ and c* H = ( m _?,l 2 _ 2H) ) 
(Here (3 denotes the Beta function.) Therefore, the map K* H defined by: 



1/2 



(K* H l [0tt] )(s) = K H (t,s)l m (s) 

is an isometry between £"(0, T) and £ 2 (0, T). This isometry can be extended to 
H(0, T), and is denoted by ^ B(<f) = ^ <j){s)dB s . 

The transfer operator K* H can be expressed in terms of fractional integrals, 
as follows: for any e £(Q,T), 

(K* H cf>)(s) = c* H T{H - l/2yi*- H I»: 1,2 {u H -Vmu)){s), 

where 

1 " T 



I%_f{s) = -f- I (u-sY^f^du 
r (a) Js 

denotes the fractional integral of / G I/ 1 (0,T), of order a G (0, 1). 

i^lf can be extended to complex-valued functions, as follows. Let £c(0,T) 
be the set of all complex linear combinations of functions l[o,t] » * G [0, r] , and 
Hc(0,T) be the closure of fic(0,T) with respect to the inner product: 



(V: ^)h c (o,t) = Qiff / / <p(tt)V'(l>)|tt-'u| 
Jo Jo 

The operator i-T^ is an isometry which maps Hc(0,T) onto L^(0, T) (the 
space of functions ip : [0, T] — » C, with J Q T | t,c ( i ) | 2 cZi < oo): for any G Wc(0, T), 

a// / T ( T <t>{u)J(v)\u-v\ 2H - 2 dudv = d H f T \l"Z 1/2 {u H - 1 / 2 ^(u))(s)\ 2 X H {ds), 
Jo Jo Jo 

(13) 

where d H = {c* H ) 2 T{H - 1/2) 2 and \ H {ds) = s l - 2H ds. 

Let £t be the class of elementary functions on [0, T) x W 1 . Note that for any 
ip G fir, the function f i-> _F<p(t, •)(£) belongs to H c (0, T), for all $ G R d . Using 
(fT2)) and (fT3j). we obtain that for any G fir, 

IM&7> = / / T I^J 172 ^- 1 / 2 ^^, •)(£))(s)| 2 A*(^)/« =: |M|g. 

JR" 1 JO 

(14) 
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We are now ready to state our result. Note that, although the conclusion 
of this result resembles that of Theorem 3 of [10] (for deterministic integrands), 
the hypothesis are different, since the proof uses techniques specific to the fBm. 

Theorem 2.1 Let [0, T] 3 t 1— > (pit, ■) £ S'{Mr) be a deterministic Junction such 
that J-tp(t, •) is a function for all t £ [0, T]. Suppose that: 

(i) the function 1 1— ► J-<p(t, •)(£) belongs to TCc(0,T) for all £ £ R d ; 

(ii) the function (t,£) i— » Tip(t, •)(£) is measurable on (0,T) x M. d ; 

(in) u"- 1 ' 2 ^ - s) h - 3 / 2 |^(m, -)(0\du < oo for all (s,£) € (0,T) x M d for 
F<p(s?-)(£) > f° r al1 ( s ,6 G ( , r ) x K<i J- 



I T :=a H / / ^(«,')K)^(v)(()h«r^«¥ffl<», (15) 

JR d Jo JO 

i/ien y £ WP arirf ||y||^-p = It- (By convention, we set (p(t, •) = /or £ > T.J 

Remark 2.2 Conditions (i)-(iii) are satisfied by the fundamental solution G 
of the wave (or heat) equation. In this case, \J-G(t, •)(£) < 1 for all i > 0, 
and hence the map t ^ !FG(t, •)(£) belongs to L^(0,T), which is included in 
«c(0,T). 

Proof: The argument is a modified version of the proof of Theorem 3.8 of [3]. 
For any £ £ R d fixed, we apply (O to the function fe(t) = F(p(t, •)(£)• We § et: 

a H / / ^(«r)(0^vKl)l«-«l aH_a d«d«= ( 16 ) 
Jo Jo 

d H i T \i%i 1/a (u B - l ' 2 r<p(u, -)(e))( s )| 2 AH(d s ). 



Jo 

It will be shown later that: 

(s,0 i-» o(«,0 := It- 1/2 (u H - 1/2 F<p(u, -)(0)(«) is measurable on (0,T) x R d . 

(17) 

Hence, we can integrate with respect to /i(d£) in (p^J|) . Using (TT5")) . we obtain: 

I T =d„ [ [ T \I?: 1/2 {u H -y 2 ^{u, ■)(0)(s)\ 2 \H(dsMd^) =: ||^||2 < oo. 
Jw d Jo 

(18) 

By the definition of HP and (Tl4|) . it suffices to show that for any £ > 0, 
there exists a function I = l £ £ £t such that: 

lk-i||o<e. (19) 

Let e > be arbitrary. By (JTTJ) and (jTHJ) , it follows that a £ L 2 ((0,T) x 
M d , Xn{ds) x (j,(d£)). Hence, there exists a simple function h(s,£) such that 

|a(*,0 - h(8,£)\ 2 \ a (d8)n(d0 < e. (20) 
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Without loss of generality, we assume that /i(s,£) — l( Ci< j](s)l^(^), with 
c,d e [0,T],c < d and A 6 S;,(M d ). By relation (8.1) of [35], we approximate 
the function l (crf] (s) in £ 2 ((0, T), \u{ds)) by l"l 1/2 {u 11 - 1 ' 2 ^))^) with Z E 
£(0,T),i.e. 

f |l(c.d](s) - I^Z ll2 {u H - l / 2 l {u)){ s )\ 2 \ H {ds) < e. (21) 
Jo 

By Lemma 3.7 of [3J, we approximate the function 1a(£) m L 2 (M. d , n(d£)) 
byThiO withZi e £(R d ), i.e. 

/ |u(0-^i(OlV^)<^ (22) 

jR d 

We define l(u,x) = l (u)li(x). Clearly I £ £ T and Tl(u, •)(£) = / (u)^ r Z 1 (^). 

Let 

b(a,£) :=/^r V ^« H - 1/2 ^(«, •)«))(*) =-frr 1/2 (« ff_1/2 io(«))(«)- ^CO- 
Usiiig d2TJ) and [|2"2]) . we obtain that: 

/ f T \h( s ,0-b( S ,0\ 2 X H (d s )fi(dO 

JR d JO 

< 2\ f / T |i (Ci<fl ( a )-j^: 1/2 (u ff - l / a Jo(«))(«)i a u(OAff(d«)/i(de)+ 

JR d JO 

< 2{e A i(A)+e[|y^ (0)r) /d i? }:=<7 1 e. (23) 
From ((20) and ([23j, it follows that 

ll¥>-i|lo = dff / / |a(s,0-6(s.OI 2 Aff(ds)/i(^) <2d ff (e + C 1 e) := C 2 e. 

JE d JO 

This concludes the proof of (fl9|) . 

We now return to the proof of (fT7|) , which uses assumptions (ii) and (iii) . If 
> 0, then (u,«,0 = li^}^- 1 / 2 ^-^^- 3 / 2 ^^,-)^) 

is measurable and non-negative, and a(s, £) = </>(u, s, £)<iii is measurable, by 
Fubini's theorem. 

Suppose next that fj u H - x / 2 {u - s) H - 3 / 2 \ Fip{u, -){C)\ du < °°. If Z(s,f) = 
l( c ,d]( s )lyl(0 is an elementary function with c,d £ [0,T],A E £>f,(R d ), then 

aiM =Irr 1/2 (« H - 1/a K«.0)(«) = U(0 f u H - l / 2 l {cA {u){u-s) H -*/ 2 du 

J s 

is clearly measurable. In general, since {u, £) i— > J-(p{u, •)(£) is measurable, there 
exists a sequence (l n ) n of simple functions such that l n (u,£) — > Tip{u 1 •)(£) for 
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all (tt,£) and \l n (u,£)\ < \Tip{u, •)(£)! f° r an ( u >£)> n ( see e -S- Theorem 13.5 of 
[5]). By the dominated convergence theorem, for every (s,£) 

|a, B («,0-o(*,OI< / u H - l ' 2 {u- S ) H -^ 2 \l n {s,O^T V {u,-m\du^Q. 



Since ai n (s, £) is measurable for every n, it follows that a(s,£) is measurable. □ 

3 The wave equation 

We consider the linear wave equation: 
d 11 

■Wp(t,x) = Au(t,x) + W(t,x), t>0,xeR d (24) 

u(0,x) = 0, xeM. d 

^(0, x) = 0, xeR d . 
dt K ' ; 

Let Gi be the fundamental solution of Uu — Aw = 0. It is known that Gi(t, •) 
is a distribution in iS'(R d ) with rapid decrease, and 

^ 1 (V)(C) = ^», (25) 
for any £ e R d , t > 0, d > 1 (see e.g. [35]). In particular, 

Gl(t,3f) = ^l{|a:|<i}, ifd=l 

Gl(M) = 2^ v / f 2_ N 2 1 ^H< t >' [id = 2 

Gi(t,x) = c d -cr t , ifd = 3, 

where <r t denotes the surface measure on the 3-dimensional sphere of radius t. 

The solution of (|2"4")1 is a square-integrable process u — {u(t,x);t > 0, x S 
R d } defined by: 

u(t,x)= / / Gi(t-s,a:-»)W(ds,dy). 

JO JR d 

By definition, u(i, x) exists if and only if the stochastic integral above is well- 
defined, i.e. g tx := G\(t — -,x — ■) G HV. In this case, E\u(t, x)\ 2 = \\gtx\\nv 
The following theorem is the main result of this article. 

Theorem 3.1 The solution u — {u(t,x);t > 0, x £ R d } of {2$ exists if and 
only if the measure [i satisfies f3|). In this case, for all p > 2 and T > 

sup sup E\u(t,x)\ p < 00, (26) 
te[o,T]xeR d 

and the map (t, x) i— > u(t, x) is continuous from R_|_ x R d into L 2 (f2). 
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Example 3.2 Let f(x) — 7a, dM ^ a ^ be the Riesz kernel of order a £ (0, d). 
Then n(d£) = |£|~ a d£ and © is equivalent to a > d - 2H - 1. 

Example 3.3 Let /(x) = 7 a / °° ^"-^/^e^e-H 2 ^ 4 ^^ be the Bessel 
kernel of order a > 0. Then = (1 + |£| 2 )~ q/2 and © is equivalent to 

a> d-2H ~1. 

Example 3.4 Let f(x) = J\%=i ( a Ht \xi\ 2Hi ~ 2 ) be the covariance function of a 
fractional Brownian field with Hi > 1/2 for all i = l,...,d. Then = 

I\Li(chM' {2H '' 1) ) and © is equivalent to £ti( 2i ^ - 1) > tf - 2H - 1. 
(This can be seen using the change of variables to the polar coordinates.) 

Remark 3.5 Condition © is equivalent to 

/ ii(d£) < oo and / U \2H+i ^ d < °°- 

J\e\<i J\i\>i \Q 2H+1 

Proof of Theorem I3.lt Note that gt x = G\(t — -,x — •) satisfies conditions 
(i)-(iii) of Theorem 12.11 Hence, g tx £ T~CP (i.e. the solution u of (f2"4")l exists) if 
and only if I t < oo for all t > 0, where 

I t :=a H / / Fg tx {u,-)(£)Fg tx {v,-)(Cj\u-v\ 2H ~ 2 dudvii{d£), 
Jm d Jo Jo 

and E\u{t,x)\ 2 = \\g tx \\ 2 nv = h- Since Tg tx {u,-){0 = e^TG^t - u, •)(£), 

JRrf JO JO 

Using (f2"5| . we obtain: 

I t= a v\ HTtP / rsin(«|?|)sin(i;|C|)| U - W | 2H - 2 d^. 

We split the integral fi(d£) into two parts, which correspond to the regions 
{|£| < 1} an< i {|£l > !}• We denote the respective integrals by Ij and l[ 2 \ 
Since the integrand is non-negative I t < oo if and only if ij: 1 ^ < oo and l[ 2 ^ < oo. 

The fact that condition © is sufficient for I t < oo follows by Proposition 
13.71 below. The necessity follows by Proposition 13.81 fusing Remark |3.5[) . 

Relation (|26p with p = 2 follows from the estimates obtained for I t — 
E\u(t, x)\ 2 , using Proposition 13.71 For arbitrary p > 2, we use the fact that 
E\u{t,x)\ p < C p {E\u(t,x)\ 2 )P/ 2 , since u(t, x) is a Gaussian random variable. 
The L 2 (fi)-continuity is proved in Proposition 13. 101 □ 

We begin with an auxiliary result. To simplify the notation, we introduce 
the following functions: for A > 0, t > 0, let 

/t(A, t) — sin rXt — r sin At, <?t(A, r) = cosrAi — cos At. (27) 
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Lemma 3.6 For any A > and t > 0, 



* 1 + A 2 ~ J R (r 2 -l) 2Uty ' ; yty ' n " * 1 + A 2 ' 
where Cj = ci(i A i 3 ) and = C2(t + i 3 ), /or some positive constants c\,C2- 
Proof: From the proof of Lemma [B. 11 we see that: 

' r [/ 2 (A,r)+. 9 2 (A,r)] = |^ ,A^)| 2 I 



(r 2 - l) 2 ' 

where <p(x) = since. Using the Plancharel's identity (|39p. we obtain: 

/ f 2 1 n2 [/ t 2 (A,r)+ g 2 (A,r)]dr= / |^ , At¥ ,(r)| 2 dr = 

f xt /"* I sin As 1 2 

27r / | sinx\ 2 dx = 27rA / | sin As| 2 d.s = 27rA 3 / - — ^ 2 ds 



The result follows using ©: (see e.g. Lemma 6.1.2) of [i3] ) 



C * 1 + A 2 - J A 2 "~ ~ ~* 1 + A 2 ' 
□ 

We denote by N t (Q the H(0, t)-norm ofun •)(£), i.e. 

N t (0 = T§f Q f sin H?D ™M£l)l« " v\ 2H ~ 2 dudv. 
Proposition 3.7 For any t > 0,£ 6 R d 



M(0 < «(t^Ip) 



1 

i 

ff+1/2 



, iei>i 



where C H = b%2 H+1 / 2 /3 and cf ] H = c h {Sh + c[ 2) )2 3ff -Va. # ere c (2) fa ^ e 
constant given by Lemma[ 



Proof: a) Suppose that |£| < 1. We use the fact that |M|^( t) — kii"IMIz,i/#(o t) — 
b 2 H t 2H ' 1 \\ip\\ 2 L2 ^ t , for any 6 L 2 (0,i), and | sina;| < x for any £ > 0. Hence, 

rt ft 



Nt{0 < b 2 H t 2H -^ I sin 2 (u\Z\)du < b 2 ^ 211 - 1 I u 2 du 

= o 2 t 2 «-lL < i 6 2 t 2ff+2 2 H+l/2 / J_ - 

ff 3 " 3 H + 
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where for the last inequality we used the fact that \ < pprep if ICI < 1- 

b) Suppose that |£| > 1. Using the change of variable u' = tt|£|, v' = v\£\, 



rt\t\ rt\6\ 

N t(0 = / / sin(u')Mv')\u'-v'\ 2H - 2 dudv 



JO 
1 ,2 



- |£|2ff+2ll sn HVim(o,t|5|)- 

Using the expression of the 7i(0, i|£|)-norm of sin(-) given by Lemma [B.ll 
we obtain: 

f |_|-(2ff-l) 

N ^) = J^T2 7 R (r2-l)2 [ftm,r)+ 9 Mlr)]dr. (28) 

We split the integral into the regions |r| < 1/2 and |r| > 1/2, and we denote 
the two integrals by iV ( (1) (£) and iV t (2) (0- 

Since / t (A,r)| < 1 + |r| and <?t(A, r) < 2 for any A > 0, r > 0, we have: 

^ } (0 < -^X,, 1/2 Tr^ [(1 + M)2 + 4]dr 



c- 



\t\<1/2 

/ -, x 2-2H -. 
CH { 1 \ 1 



1-H\2J |£| 2ff+1 ' 

We used the fact that |£| 2ff + 2 > \£\ 2H+1 if |f| > 1, and [(l + |r|) 2 +4] < 

I 3^ F [(3/2) 2 +4] = Cif |r|<l/2. 

Using the fact that | T |-( 2 ^-i) < (i)-( 2 ^-i) if | T | > | ; LemmaEH and the 
fact that |£| 2 /(i + |£| 2 ) < 1, we obtain: 

* 2 -g-D ^ / T| , 1/2 (^TF [/t2(iar)+gt2(lel ' r)]rfT 

< cg c ( 2 ) 1 |£| l^' 2 
- 2 -( 2 ^-!) ' |£| 2ff + 2 ISI 1 + |£| 2 

< Cg J 2 ) 1 



2~(2H-1) "t |£|2.ff+l' 
□ 

Proposition 3.8 a) If It < 00 / or * = 1; fieri J|£|<i M^O < °°- 

6) Lei I > 1 be the integer from 110\) and m = 21 — 2. For an?/ i > 0, 



/" Kg) ^ „ /v- »n r (2) , m +i /" , 9Q x 
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where ajj,t = % 2H J '(chc[ ), b t — 2C/q and c\ is the constant of Lemma \3.6\ 
In particular, if ij® < oo for some t > 0, then fi^i CI ^^ 2H+1 ^ fJ-(d^) < oo. 

Proof: a) Using the fact that sinx/x > sin 1 for all x £ [0, 1], we have: 

I« = f f (\in{u\Z\)sm{vmu-v\™-*dudv 

J\i\<i Is I Jo Jo 

> sin 2 l f fi(d() [ [ uv\u-v\ 2H - 2 dudv. 
J\Z\<i Jo Jo 

b) According to ([25]). 

4 2) = «* J 0& J ^™~y [m\,T)+gX\Z\,T)]dT. (30) 
For any fc G { — 1, 0, . . . , m}, let 

By (nnn, /(m) = s mi ier (2H+2+m V(^) < 4^ i^r 2 V(^o < 00. 

We will prove that the integrals /(fc) satisfy a certain recursive relation. By 
reverse induction, this will imply that all integrals /(fc) with fc € {—1, 0, . . . , m} 
are finite. For this, for k € {0, 1 . . . , m}, we let 

Mk) :=^ >i] J^ | KF l_^[ /t2( |e|,r)+ 5 2 (|C|,r)]dr. (31) 

We consider separately the regions {|r| < 2} and {|r| > 2}. For the region 
{|r| < 2}, we use the expression ([3D} of / t (2) . Using the fact that |£|2ff+2+fc > 
|C| 2H+2 (since fc > 0), and \ T \-VH-i) > 2-( 2ff - 1 ) if | T | < 2, we obtain: 

m) - / Elil KF^l S2 F 1 TF [/ ' !(laT) + 9 ' <H -* T 

< 2 2H - 1 — /f, by®. 

For the region {|r| > 2}, we use the fact |/t(A, r)| < 1 + r and |<?t(A, r)| < 2 
for all A > 0, r > 0. Hence, 

4'M : = / £l>1 i«F^/ T|>! F ?1 TF l/?(l?l,T) + ! ' f(l?l,T)l ' iT 
s L^L^V 1 + H)2 + 4| * = c/(t) ' 
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Hence, for any k € {0, 1, . . . , to} 

A t (k) < 2 2H - 1 —4 2) +CI(k). 

CH 

Using Lemma [3761 and the fact that jjrigp > § if |£| > 1, we obtain: 

for all A; G {0, 1, ... , to}. From the last two relations, we conclude that: 

U 1) I(k-l)<2 2H - 1 — 4 2) +CI(k), Vfce{0,l,...,m}, (32) 
2 c H 

(2) 

or equivalently, /(fc — 1) < a,H,t!t + btl(k), for all A; € {0, 1, . . . , to}. Relation 
(j2"9")l follows by recursion. □ 

Remark 3.9 In the previous argument, the recursion relation (|32[) uses the fact 
that k is non-negative (see the estimate of A' t {k)). Therefore, the "last" index 
k for which this relation remains true (counting downwards from to) is k = 0, 

leading us to the conclusion that J^ >1 |£| ^ H+1 V(^£) < °°> if H < 00 • 

The next result shows that the map (t, x) — > u(i, x) from R+ x M d into L 2 (f2) 
is continuous. 

Proposition 3.10 Suppose that $8}) holds, and let u = {u(t,x),t >0,xG R d } 
be the solution of {21$ . For any t > 0, 

E\u(t + h, x) — u(t, x)\ 2 — > as \h\ — > 0, uniformly in x G R d (33) 

and 

- u(i,y)| 2 -> as \x — y\ — » 0. (34) 

Proof: We use the same argument as in Lemma 19 of [10 (see also the erratum 
to [10]). We first show (|33|. 

Suppose that h > 0. Splitting the interval [0, t + h] into the intervals [0, t] 
and [t, t + h], and using the inequality \a + b\ 2 < 2(a 2 + b 2 ), we obtain: 

E\u(t + h,x) -u(t,x)\ 2 < 2{\\(g t+ h,x - gtx)l[o,t\\\nv + \\9t+h,xl[t,t+h] \\nv) 

=: 2[E 1>t (h) +E 2 (h)]. 



Since H9t+h, x -9tx){u, 0(0 = e-* x FGi(t + h-u, 0(0 - .FGi(t - u, 0(0, 



t ct 



E ltt (h) = a H Kd£) / dvdv\u-v\ 21i - 2 T(g t+h , x -g tx )(u r )(0 



f{gt+h,x- gtx)(v, 0(0 

oh / **(#)/ / dudv|u-tf ff ~ 2 [.FGi(u + /v)(0 -JR7i(tt,0(O] 



^d(i; + /i,0 (0-^(^,0(0 
/x(dO 



i^| 2 -*t(^lfl 



15 



where 
kt{h,\£\) 



(sin((u + h)\£\) - sin(u|C|))(sin((« + - sin(«|S|)) 



JO 
|2/f-2 



|u - v\ 2 "- 2 dudv = || sin((- + h)\i\) - sin(- |£|)|& (0)t) . 



By the Bounded Convergence Theorem, lim^ |£|) = 0, for any £ e M d . 

The fact that Ei_ t (h) — > as /i J, will follow from the Dominated Conver- 
gence Theorem, once we prove that: 



h(h, Ifl) < fctdfl), V/ie [0,l],V£e 



and 



lei 2 



fe t (|£|) < co. (35) 



When |£| < 1, using the same argument as in Proposition 13.71 we get: 
k t (h,\Z\) < b^-'Wsmd- +h)\Z\)-M- m\l H o,t) 



sin {{u + h)\£\)du + / sm\u\£\)du 



{u 2 + l)du+ / u 2 du) =:k t (\£\) 



< 2b 2 H t 2H - 1 \^\ 2 ^2(u 2 + l)du + J 

Suppose that |£| > 1. We use the fact that: 

k t (h, |e|) < 2(|| sin((- + h)\Z\)\\n m + II M- l£l)ll«(o,t))- 

Using the change of variables u' = (u + h)\£\,v' = (v + h)\£\, and ([40 
(Appendix A) we obtain: 



l|sin((- +h)\Z\)\\nm 



OtH 

OiH 
\i? H 

C H 
2H 



sin((u + /i)|£|) sm((v + h)\t\) \u - v\ 2H ~ 2 dudv 

(t+h)\t\ r(t+hM\ 
h\£\ Jh\(\ 



sin(u') sin(v')|u' - v'\ 2H ~ 2 du'dv' 



lei 



where tp(t) — suit. Note that the square of the real part of J-'} l \£\(t+f l -\i£iip(T) is: 



(t+h)\£\ 



cos rt sin tdt 



2 




< 2 


L 







<t+h)\s\ 



cos rt sin tdt 



cos rt sin tdt 



and the square of the imaginary part of ? r h\£\,(t+h)\£\ l P{' r ) i s: 



r(t+fc)|C| 



sin rt sin tdt 



2 




2 


< 2 


/ sin rt sin idi 


+ 2 




Jo 





sin rt sin tdt 
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We now use the following fact (see Appendix B): for any T > 



cos rt sin tdt 



sin rt sin tdt 



1 



(r 2 - l) 2 

From here, it follows that k t (h, |£|) is bounded by 
2c H 



[(sin tT-t sin T) 2 +(cos rT-cos T) 2 



2H 



-(2H-1) 



(r 2 - 1) 



^[/t+/ l (lC| J r)+9 t V(lei,T)+^(|^,T)+^(|^,T)+//(|e|,r)+^(|C|,r)]dr 5 



where /t(A, r) and gt(A,r) are defined by 1|27|) . The argument of Proposition 
13.71 shows that for any i > and |£| > 1 



| r |-(2ff-l 



(4) 



•; 2 _ 1)2 [^(l^l,r)+g?(|^|,r)]dr < C^j^ja ^ 4?M 

where c| 4 ^ = (I) 2 2H _|_ f |) ^ 2H ^ c | 2 \ Since c[ 4 ^ is non-decreasing in t 
and /i S [0, 1], fc t (/i, |£|) is bounded by 

2Cg (- p W r W I r W N < 2c g f „(4) (4) r (4)x._ fr /|> n 

|£|2fl-l ^t+hM i c h,H ' c t,H> — |£|2ff-l V L t+l,if T C 1H f c ( H j .— ftt(J?|,/- 



This concludes the proof of (|35|) . 

A similar argument shows that E2Q1) — ► as h [ 0, since 



£ 2 (/i) 



t+/i r t+h 



TG\(t + h-u, ■){^)J : G l {t + h-v, -)(0I« - u| 2ff ~ 2 dudw/i(dO 

I?! 2 



Jo 



sin(u|£|) sin(t>|f - v\ 2H ^ 2 dudv 



The case h < is treated similarly. Using the same argument as above, it 
follows that for any h > 0, E\u(t - h,x) - u(t, x)\ 2 < 2{E' l t {h) + E 2 {h)), where 

KtW = I ^$-K(h, |£|), and k' t (h, = || sin(- |£|) -sin((. - fc)l£l)llW 

To prove (|34j) . note that 
£|u(i, sc) - u(t, y)| 2 = ||fe - gtyWnv = 



a H 



Jo 



F{9tx ~ 9t y ){u, -)(C)^(9tx ~ 9t y ){v, -){£,)\u - v\ dudv/j,(d£,) 

•t r t 



Jo 



FGx{u, -){i)TGi{v, ■)(£)!« - v \ 2H ~ 2 dudvn(dO, 



which converges to as \x — y\ — > 0, by the Dominated Convergence Theorem. 
□ 
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Example 3.11 There exists an interesting connection between the solution of 
the wave equation with fractional noise in time and Riesz covariancc in space 
and the odd and even parts of the fBm. Indeed, if / be the Riesz kernel of order 
a E (0, d), then 

It = a H [ dm-"- 2 "- 2 I TT^y [f!m,r)+g!m,r)}dr 

f \t\- {2H - 1] ( f°° (smT\t-Tsm\t)\ g „ f°° (cosrAi-cos 
= 2a H c d / — T^r / " ^ -\~ e d\ 1 ' 



(r 2 I) 2 \Jo y 

where 6* = a + 1 — d + 2 if > under ([5]). If 8 < 1, the two integrals e?A can be 
expressed in terms of the covariance functions of the odd and even parts of the 
fBm (see QU). 



4 The heat equation 

In this section, we consider the the heat equation with additive noise: 

^(*,a0 = ^Au(t,x) + W(t,x), t>0,xeR d (36) 
u(0,x) = 0, x E M d . 

Equation (|3^|) was treated in (3], in the case of particular covariance kernels 
/. We give here an unitary approach which covers the case of any covariance 
kernel /, which satisfies ([9])- 

The case of the heat equation is actually much simpler than the case of the 
wave equation, since both the fundamental solution G and its Fourier transform 
are non- negative functions. 

More precisely, let Gi be the fundamental solution of u t — \ Au = 0. Then 

G 2 (M) = ^cx P (-^), ;>o,,e^ 

and 2 

jPG 2 (V)(0 = expf-^|Lj, t>0^ER d . (37) 

We will prove the following result. 

Theorem 4.1 The solution u = {u(t, x),t > 0, x E R d } of iS6}) exists if and 
only if the measure fi satisfies {P^. In this case, i26\) holds for all p > 2 and 
T > 0, and the solution is L 2 (fl)- continuous. 

Remark 4.2 (i) When / is the Riesz kernel of order a, or the Bessel kernel of 
order a, condition © is equivalent to a > d — AH . When / is the covariance 
function of the fractional Brownian field with Hi > 1/2 for all i = l,...,d, 
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condition ^ is equivalent to Yli=i(2Hi ~ 1) > d — AH. Note that this condition 
is weaker than the condition given in [33) . 

(ii) In Theorem 2.1 of the Erratum to [3] it has been proven that condition 
© implies that ||flta;||w < 00 I0r an y t >0 and x G M. d . 

Proof of Theorem I4.lt Note that g tx = £?2(f — -,x — •) is non-negative. 
Hence, gtx £ HV if and only if g tx e \TCP\. This is equivalent to saying that 
Jt \\9tx\\fn-p\ < 00 for all t > 0. Note that 

Jt = a H / g tx {u,y)gtx{v,z)f(y - z)\u- v\ 2H ~ 2 dydzdudv 

Jo Jo J& d JR d 
t ft 



a H I ^9tx{u,-){S)Tg tx {v,-){C)\u ~ v\ [i{d£)dudv 

Jo Jo 



t ft 



a H TG 2 {t-u,-){i)TG2{t-v,-){i)\u~v\ IH -^{di)dudv. 



o Jo 

Using (|37|) and Fubini's theorem (whose application is justified since the 
integrand is non-negative), we obtain: 

J t = an J J* J* exp (~^f-) exp (-^) I" - v\ 2H - 2 dudv 

The existence of the solution follows from Proposition 14.31 below, which also 
gives estimates for J t = E\u(t, x)\ 2 (and hence for E\u(t, x)\ p ). The L 2 (Q)- 
continuity is given by Proposition 14.41 □ 

Let 

2\ / „,l£|2\ 

2H-2, 



\u-v\ ZH ~ z dudv 



The next result is similar to Lemma 6.1.1) of [43 . 
Proposition 4.3 For any t > 0, £ G M d , 

where C H = b 2 H (AH) 2H . 

Proof: Suppose that |£| < 1. Using the fact that IMIw(ot) — fr/f* 2if_1 IMIz,2(o t) 
for all <p e L 2 (0,t), e~ x < 1 for any x > 0, and \ < if |£| < 1, 

1 

~\ 



MO < b 2 H t 2H ~i J* eM-u\e)du < b 2 H t 2H < b 2 H 2 2H t 2H 
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Suppose that |£| > 1. Using the fact that ||y||^( t) — kj?IMIi,i/ff(o t) ^ or an y 
<p e L^ H (0,t), 1 - e~ x < 1 for all x > 0, and < i+rcp i we obtain: 



exp 



^t(0 < & 

< 6|r(4H) aH 



2# 
1 



2ff 



2ff 



2ff 
IP 



2ff 



1 — exp 



2H 



2H 



This proves the upper bound. 

Next, we show the lower bound. Suppose that i|£| 2 < 1. For any u £ [0, t], 
Hl|L < ML < I Using the fact that e~ x > 1 — x for all x > 0, we conclude 
that: 



exp 



> . "lei 2 i w c rn 

> 1 — ^— > 2' Vtt e [M- 



Hence 

A t (0 > a H 



2 „t r i 



JO 



\u - v\ 2H - 2 dudv = li 2H > V H 
11 4"4 



2ff 



For the last inequality, we used the fact that 1 > 



Suppose that i|£| 2 > 1. Using the change of variables u' = u|£| 2 /2,?/ 



v\(\ /2, we obtain: 



2 2« /■ 



Since the integrand is non-negative, 

22H rl/2 

A t {0 > a Hj7m 



9 2H r m 2 /2 r t\a\ 2 /2 

1 e- u e- v \u'-v'\ 2H - 2 du'dv'. 



lei 4 



1 



e- u e- v \u-v\ 2H - 2 dudv 

2H+2 



H(0,l/2) jljiH 



> 2 



2ff 



2ff 



where for the last inequality we used the fact that ^ > pprcra , and 



(|) 2ff+2 . (This follows since e" 11 > 1 - u > \ for all u G [0, |]. 



l«(0,l/2) 



> 



□ 



Proposition 4.4 Suppose that (0|) holds, and let u = {u(t,x),t > 0, x 6 M d } 
6e i/ie mild-sense solution of [3b]) . Then the map (t,x) — > u(t,x) from K+ x M d 
mio L 2 (f2) zs continuous. 

Proof: The argument is similar to that of Proposition 13.101 In this case, if 
h> 0, 

£i,t(fc)= / M(^)fc t (^, 1^1), 
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where 




2 dudv. 



We omit the details. □ 



Remark 4.5 We consider the operator Lu = dtu—J2i j—i a ijdx i x j u ~J2i=i bid Xi u. 
Let Gs(t, x; s, y) be the fundamental solution of Lu = 0. We assume that: 

(i) The functions a tJ ,b t : [0, T] x R d -> R, i, j = 1, . . . , d are a/2-H61der 
continuous in i and a-H61der continuous in x, for some a € (0, 1). 

(ii) There exist some fc, X > such that for all (t, x) E [0, T] xI^ie R d , 

d 

m 2 < J2 a<At,x)tiZi<K\Z\ 2 - 

Under these assumptions, G3 is a positive function defined on [0, T] x R d x 
[0, T] x R d n {(s, i); < s < t < T}, which satisfies: (see p. 376 of [21]) 

G 3 (t,x;s,y) < c x {t - s)~ d ' 2 exp (-c 2 ^~^ j := G' 2 (i - s, a; - y). (38) 

Since G' 2 (t, x) is essentially the same as the heat kernel G2(t, x), the solution 
of Lu(t, x) — W(t, x) (with vanishing initial conditions) exists, if the measure \i 
satisfies condition ((9]). 



A Some useful identities 

Recall that the Fourier transform of a function ip 6 L X (R) is defined by: 

Tip{r) = [ e- iTX (p(x)dx. 
Jm 

For an interval (a, b) C R, we define the restricted Fourier transform of a function 
Lp S L 1 (a.b): 




(p(x)dx = T{<pl [aM ){T). 



One can prove that Ftp £ L 2 (R), for any tp e L 1 (R) n L 2 (R). By the Plan- 
charel's identity, for any tp,ip g L 1 (R) n L 2 (R), we have: 

/ p(x)^(a;)da; = (27T)" 1 / ^(t)^(t)^. 
7r jr 
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In particular, for any ip, ip £ L 2 (a, b), we have: 



b 



<p(x)ij;(x)dx = (27T)- 1 / FaM T )F*,bMT)d£. (39) 



(Consider <p = tpl [a ,b]- Then <p e L 1 (R) n L 2 (R) and = F a ,b<p(0-) 

The proof of Theorem l3.1l uses. in an essential way, a formula for the H(0, T)- 
norm of sin (developed in Appendix B), which is in turn based on the following 
result. (This result can be derived using for instance, the results of [34].) 

Lemma A. 1 Let H E (± 1). For any ip,ip € L 1 (R) n L 2 (R), 



Lp{u)^{v)\u-v\ 2H - 2 dudv = c H / T(p{T)Ti}{T)\T\-^ H -^dT 



where a H = H(2H - 1) and c H = T(2H + 1) 8m(irH)/(2ir). 
In particular, for any ip,ip E L (a,b), 



b ,-b 

<p(u)^(v)\u-v\ 2H - 2 dudv = c H I F a ,b<p(r)F a ,biKT)\T\-l 2H - 1 '>dr. 



(40) 



B The H(0, T)-norm of sin 

Lemma B.l Let tp(t) = suit, t E [0,T]. Then 

.. a f (sin tT-t sin T) 2 + (cos tT - cosT) 2 ,_ (2 ir_i). 

IMIh(0,T) = C * / R - p _ 1)2 M ( W > 

where c H = T(2H + 1) sm(TrH) / (2w) . 
Proof: By gDJ), 

2ff-2, ,„_^„ /"|^_, / x,2, |-(2H-l: 



IMIwfo.r) = / / <p(«Mv)|w-u| dudv = c H / IJ^o.t^Ct)! |r| 

JO JO JR 

Note that |-Fo,!t<p(t)| 2 = J Q T e- lTt <p(t)dt * = I 2 + J 2 , where 

,T ,T 

Ji = Rc^o.t^t)] = / cos Tt sin tdt, J\ = Im[JFo,T</?(T)] = I saxTtsmtdt. 
Jo Jo 

We calculate I\ first. Using integration by parts, we obtain: 

h = 1 -costTcosT -tI 2 , 
where I2 = J Q T sin rt cos tdt. On the other hand, 



h + h= [ sin[(r + 
Jo 



1 - cos[(r + 1)T] 



r+1 
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Solving for 7i and I2, we obtain 
1 

1 -T 

Similarly, letting J2 = J Q T cos rt cos tdt, obtain: 



h = — (1 — cos tT cos T — t sin tT sin T) . 



r J 2 - Ji = sin rT cos T and J 2 - J x = — !— sin[(r + l)Tl. 

T + 1 



Solving for Ji, we obtain: 



Ji = - 1 , (r cos rT sin T - sin tT cos T) . 

1 — T 2 



An elementary calculation shows that: 

7 1 + J l = M 1 K Sm rT - T Sm T f + ( COS TT - C0S T ) 2 ] ' 

(1 — t^y 

□ 

Remark B.2 Let B = (B t )teR a fBm of index H (on the whole real line). Let 
B° = {B°) tm and B e = {B^) tm be the odd and even parts of B (see [IS]). B° 
and _B e are independent centered Gaussian processes with B t = B° + Bf , and 

E{B°B° S ) = (l m ,l {0iS) ) o :=c H [ ^^ ^-m-l)^ 

Jr t 

p/™ Rel ,, 1 \ / (1 - cos rt)(l- costs) (2H _ 1} 

h \ B t B s) - \ 1 (o,t)' 1 (o.s))e ■- c H J — 2 \t\ 'dr. 

In general, for 99 e i 1 (M) n L 2 (R), we have: 

£[B°^) 2 ] = IMIo := <=h / \R C [T V (tW\t\-^ h - 1 Ut 

Jr 

25[B» 2 ] = Ml := c ff / llm^MlHrr^- 1 )^ 

Jr 

In the proof of Lemma IB. II 1\ = Ii (r) and J\ = J% (r) are the real and 
imaginary parts of f(iplt w)(r), where ip(t) = sini. Hence 

E 



/ smtdB?) =|| S m(-)l [ o,T]||^ = c H [ \h(r)\ 2 \r\-^ H -^dr 
Jo J Jr 

E(£wLtdBz\ =|| S in(.)l [ o,T]||^ = c H JjMT)\ 2 \r\-^ H -^dT. 



Acknowledgement. The authors would like to thank Professor Robert Dalang for the 
invitations to visit EPFL and the useful discussions. 



23 



References 



[1] Alos, E., Mazet, O. and Nualart, D. (2001). Stochastic calculus with respect 
to Gaussian processes. Ann. Probab. 29, 766-801. 

[2] Alos, E. and Nualart, D. (2003). Stochastic integration with respect to the 
fractional Brownian motion, Stoch. Stock. Rep. 75, 129-152. 

[3] Balan, R.M. and Tudor, C. A. (2008). The stochastic heat equation with 
fractional-colored noise: existence of the solution. Latin Amer. J. Probab. 
Math. Stat. 4, 57-87. Erratum in Latin Amer. J. Probab. Math. Stat. (2009) 
6, 343-347. 

[4] Balan, R.M. and Tudor, C. A. (2009). Stochastic heat equation with multi- 
plicative fractional-colored noise. To appear in J. Theoret. Probab. Preprint 
available on larXiv:0812.1913[ 

[5] Billingsley, P. (1995). Probability and Measure. Third Edition. John Wiley, 
New York. 

[6] Biagini, F., Hu, Y., 0ksendal, B. and Zhang, T. (2008). Stochastic Calculus 
for the Fractional Brownian Motion and Applications. Springer. London. 

[7] Carmona, R. and Nualart, D. (1988). Random nonlinear wave equations: 
Smoothness of the solutions. Probab. Th. Rel. Fields 79, 469-508. 

[8] Conus, D. and Dalang, R. C. (2009) The non-linear stochastic wave equa- 
tion in high dimensions. Electr. J. Probab. 22, 629-670. 

[9] Dalang, R. C. and Hou, Q. (1997). On Markov properties of Levy waves in 
two dimensions. Stoch. Proc. Appl. 72, 265-287. 

[10] Dalang, R. C. (1999). Extending martingale measure stochastic integral 
with applications to spatially homogenous s.p.d.e.'s. Electr. J. Probab. 4, 
no. 6, 29 pp. Erratum in Electr. J. Probab. 6 (2001), 5 pp. 

[11] Dalang, R. C. and Frangos, N. E. (1998). The stochastic wave equation in 
two spatial dimensions. Ann. Probab. 26, 187-212. 

[12] Dalang, R. C. and Sanz-Sole, M. (2008). H older- Sobolev regularity of the 
solution to the stochastic wave equation in dimension 3. Memoirs of the 
AMS. 

[13] Da Prato, G. and Zabczyk, J. (1992) Stochastic Equations in Infinite Di- 
mensions, Cambridge University Press. 

[14] Decreusefond, L. and Ustunel, A. S. (1999). Stochastic analysis of the frac- 
tional Brownian motion. Potential analysis. 10, 177-214. 

[15] Dzhaparidzc, K. and van Zanten, H. (2004). A series expansion of fractional 
Brownian motion. Probability Theory and Related Fields. 130, 39-55. 



24 



[16] Gradinaru, M., Russo, F. and Vallois, P. (2003). Generalized covariations, 
local time and Stratonovich It's formula for fractional Brownian motion 
with Hurst index H>\. Ann. Probab. 31, 1772-1820. 

[17] Gubinelli, M., Lejay, A. and Tindcl, S. (2006). Young integrals and SPDEs, 
Potential Anal. 25 (2006) 307-326. 

[18] Hu, Y. (2001). Heat equations with fractional white noise potentials. Appl. 
Math. Optim. 43, 221-243. 

[19] Hu, Y. and Nualart, D. (2009). Stochastic heat equation driven by frac- 
tional noise and local time. Probab. Theory Rel. Fields 143, 285-328. 

[20] Krylov, N. V. (1999). An analytic approach to SPDEs, in: Stochastic partial 
differential equations: six perspectives, Math. Surveys Monogr. 64, 185-242, 
AMS, Providence, PJ. 

[21] Ladyzhenskaia, O. A., Solonnikov, V. A. and Ural'tceva, N. N. (1968). Lin- 
ear and Quasi-linear Equations of Parabolic Type. Translations of Mathe- 
matical Monographs 23, AMS, Providence. 

[22] Lyons, T. (1998). Differential equations driven by rough signals, Rev. Mat. 
Iberoamericana 14, 215-310. 

[23] Lyons, T. and Qian, Z. (2002). System Control and Rough Paths, Oxford 
University Press. 

[24] Maslovski, B. and Nualart, D. (2003). Evolution equations driven by a 
fractional Brownian motion. J. Fund. Anal. 202, 277-305. 

[25] Millet, A. and Sanz-Sole, M. (1999). A stochastic wave equation in two 
space dimensions: Smoothness of the law. Ann. Probab. 27, 803-844. 

[26] Mishura, Y. S. (2008). Stochastic Calculus for Fractional Brownian Motion 
and Related Processes. Springer. Berlin. 

[27] Mueller, C. (1997). Long time existence for the wave equation with a noise 
term. Ann. Probab. 25, 133-151. 

[28] Nualart, D. (2003). Stochastic integration with respect to fractional Brow- 
nian motion and applications. Contemp. Math. 336, 3-39. 

[29] Nualart, D. (2006). The Malliavin Calculus and Related Topics. Second 
Edition. Springer. Berlin. 

[30] Nualart, D. and Quer-Sardanyons, L. (2007). Existence and smoothness 
of the density for spatially homogenous SPDEs. Potential Anal. 27 (2007) 
281-299. 

[31] Nualart, D. and Vuillermont, P.-A. (2006). Variational solutions for partial 
differential equations driven by fractional a noise. J. Fund. Anal. 232 390- 
454. 



25 



[32] Nualart, E. and Viens, F. (2009) The fractional stochastic heat equation 
on the circle: Time regularity and potential theory. Stock. Proc. Appl. 119, 
1505-1540. 

[33] 0ksendal, B. and Zhang, T. (2001). Multiparameter fractional Brown- 
ian motion and quasi- linear stochastic partial differential equations. Stock. 
Stock. Rep. 71, 141-163. 

[34] Pipiras, V. and Taqqu, M. (2000). Integrations questions related to frac- 
tional Brownian motion. Probab. Theory Rel. Fields 118, 251-291. 

[35] Pipiras, V. and Taqqu, M. (2001). Are classes of deterimistic integrands for 
the fractional Brownian motion on a finite interval complete? Bernoulli 7, 
873-897. 

[36] Quer-Sardanyons, L. and Sanz-Sole, M. (2004). Absolute continuity of the 
law of the solution of the 3-dimensional stochastic wave equation. J. Fund. 
Anal. 206, 1-32. 

[37] Quer-Sardanyons, L. and Sanz-Sole, M. (2004). A stochastic wave equation 
in dimension 3: Smoothness of the law. Bernoulli 10, 165-186. 

[38] Quer-Sardanyons, L. and Tindel, S. (2007). The 1-d stochastic wave equa- 
tion driven by a fractional Brownian sheet. Stock. Proc. Appl. 117, 1448- 
1472. 

[39] Sanz-Sole, M. and Sarra, M. (2000). Path properties of a class of Gaussian 
processes with applications to spde's. Canadian Math. Soc. Conf. Proc. 28. 

[40] Sanz-Sole, M. and Sarra, M. (2002). Holder continuity for the stochastic 
heat equation with spatially correlated noise. In: "Progess in Probab." 52, 
259-268, Birkhauser. 

[41] Saint-Loubert Bie, E. (1998). Etude d'une EDPS conduite par un bruit 
poissonnien. (Study of a SPDE driven by a Poisson noise). Probab. Th. 
Rel. Fields 111, 287-321. 

[42] Sanz-Sole, M. and Vuillermont, P.-A. (2009). Mild solutions for a class of 
fractional SPDE's and their sample paths. J. Evolution Equations 9, 235- 
265. 

[43] Sanz-Sole, M. (2005). Malliavin Calculus with Applications to Stochastic 
Partial Differential Equations. EPFL Press, Lausanne. 

[44] Tindel, S., Tudor, C.A. and Viens, F. (2003). Stochastic evolution equations 
with fractional Brownian motion. Probab. Th. Rel. Fields 127, 186-204. 

[45] Treves, F. (1975). Basic Linear Partial Differential Equations. Academic 
Press, New York. 



26 



[46] Walsh, J. B. (1986). An introduction to stochastic partial differential equa- 
tions, Ecole d'Ete de Probabilites de Saint-Flour XIV, Lecture Notes Math. 
1180, 265-439, Springer, Berlin. 

[47] Zahle, M. (1998). Integration with respect to fractal functions and stochas- 
tic calculus I, Probab. Th. Rel. Fields 111, 333-374. 



27 



